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Single-instanton approximation (SIA) is often used to evaluate analytically instanton contributions 
euclidean correlation function in QCD at small distances. We discuss how this approximation can be 
consistently derived from the theory of instanton ensemble and give precise definitions to a number 
of different "quark effective masses" , generalizing the parameter m* , which was introduced long ago 
to account for the collective contribution of the whole ensemble. We test numerically the range of 
applicability of the SIA for different quantities. Furthermore, we determine all the effective masses 
(for random and interacting instanton liquid models) as well as from phenomenology, and discuss 
to what extent those are universal.lan 
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I. INTRODUCTION 



The instanton liquid model of the QCD vacuum |^ is based on a semiclassical approximation, in which all gauge 
configurations are replaced by an ensemble of topologically non-trivial fields, instantons and anti-instantons. It 
T-H " remains a model because we do not still understand why large-size instantons are not present in the ensemble. Fits to 
CO phenomenology and later lattice studies showed that their total density is no ~ 1/to~* while the typical size of about 
p ^ 1/3/to, leading to small dilutcness parameter nop'^ ^ 10^^ With these parameters, the model quantitatively 
^ . explains such important phenomena as spontaneous SU{Nf) chiral symmetry breaking for Nf quark flavors, the 
Qv^ ■ explicit U(l) symmetry breaking, and many more other details of hadronic correlators and spectroscopy (for a recent 
, example see discussion of vector and axial correlators for a review see [^). The main feature of the instanton 

■ ensemble is that each pseudo-particle is an effective vertex with 2Nf quark lines js), which are exchanged between 
', them and fill the vacuum. A theory is developed, called Interacting Instanton Liquid Model (IILM) which include 
' these 't Hooft interactions to all orders Q. 

\^ '\ If new sources (external currents) are added, they produce extra quarks which interact with those in vacuum and 

■ produce non-trivial correlation functions. In particular, many (Lorentz scalar) chirally odd local operators obtain 
' non-zero vacuum expectation values. In general, all of those "condensates" and correlation functions are determined 

by the interaction of instantons and thus depend on the global (collective) properties of the ensemble. 
^ ' On the other hand, as the instanton vacuum is fairly dilute, one may think that the correlation functions at distances 
D , short compared to instanton spacing x R — n~^/'^ ~ l/^ri may be dominated by a single instanton, the closest 
i-S^ (or leading) one (LI). This framework ( which we shall refer to as the "single instanton approximation", SIA ) has 
the advantage to allow to carry out calculations analytically. It is therefore possible to obtain closed expressions for 
instanton contribution to Green's functions in momentum or in Borel space. 

In SIA collective contribution of all instantons other than the leading one are taken care of by a single effective 
parameter, usually called effective mass, m*. In the simplest approximation, it can be associated with an average 
value of the quark condesate [Q : 

2 

m* — m TT^ p^ < uu >, (1) 

3 

which leads to the value m* ~ llOMeV Note that it is already very different from what one infers from the same 
model for long distance (or zero Euclidean momentum) limit of the quark propagator, which gives constituent quark 
mass of the order of 400 MeV. 

Furthermore, although the SIA has been used in several phenomenological studies (e.g. [|], ^, 0, and references 
therein), its derivation was never discussed in detail, its range of applicability was never quantitatively checked, and 
the values of relevant effective masses well specified. And indeed, if one uses the value m* ~ llOMeV the correlation 
functions, evaluated in the SIA, do not agree with the results of the random and interacting instanton liquid 

In this paper we identify the origin of such discrepancy and calculate the values of effective mass appropriate for 
different observables. This analysis reveals that the discrepancy between SIA and full liquid calculations is due to 
an incorrect estimate of the effective mass, to*. We also present a systematic study of the SIA in QCD by itself. 
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We show that the approach is really accurate only for calculations that involve operators of dimension six or more, 
or correlators with more than one zero-mode propagator. We shall also prove that the mass terms, appearing in 
matrix elements involving different numbers of zero-mode propagators, are indeed independent parameters that have 
to be fixed separately. We provide with the definitions of all such mass factors in terms of averages of the instanton 
ensembles and prove that they are nearly universal, i.e. the same for all similar correlation functions. 

The paper is organized as follows. In section |ll| we derive the SIA from the theory of the instanton ensemble, 
in section III we present the results of our numerical simulations that estimate the contribution from the leading- 



instanton to several correlation functions. In section tV, we evaluate the effective mass terms both from the random 



and interacting instanton liquid and compare it with the values obtained phenomenologically from the pion sum-rule. 
In section ^ we compare our effective masses with the so-called "determinantal masses" , which are other effective 
parameters that can be defined in terms of averages of the fermionic determinant. The main results of our analysis 
are summarized in section VI. 



II. QUARK PROPAGATOR 



In this section we review how the quark propagator in the instanton vacuum is obtained and present consistent 
derivation of the SIA. 

The quark propagator in general background field is 



Si{x,y) =< x\{il^j + im) >, 



(2) 



where I/>i denotes the Dirac operator. The inverse (|^) can be formally represented as an expansion in eigenmodes of 
the Dirac operator: 



Siix,y) = 



X + im ' 



ilpiipx{x) = Xip\{x). 



(3) 



From eq. it follows that the propagator of light quarks is dominated by eigenmodes with small virtuality. 

We begin by considering the academic case in which the vacuum contains only one isolated instanton. One eigenmode 
of with zero virtuality (zero-modes) is given by 't Hooft |^ , : 



ii?Mx) = 0, 



p 



■tpOaA^i z) = - 



1 



TT ((x - z)2 + p2)3/2 



1-75 



2 vO^^ 



(4) 



Q/3 



where z denotes the instanton position, a, (3 = 1, • • - 4 are spinor indices and Uat represents a general group element. 
Isolating the contribution from zero-modes we can write: 

a, X - z) i^liy - z) ^ tpxix- z)^Pl{y~ z) 

^nXiViZ) — : h > — — 

t m — ' A + tm 

= Sr{x,y;z) + Sr"'{x,y-z). (5) 
The zero-mode part of the propagator in the field of one instanton can be evaluated from (||) and to give, ||l0|] : 



5r(x,y;z) 



8m 



1-75 



{x- z)(l3{y- z), 



where 



1 



7r\t\{t^+p^) 

The corresponding expression in the field of one anti-instanton is obtained through the substitution: 

1 - 75 1 + 75 - + 



(6) 



(7) 



(8) 



3 



In the chiral limit, m 0, the expression for S'"^'"(x, y; z) is also known exactly [jll]. In the limit of small distances 
(|a; — 2/1 — > 0), or if the instanton is very far away (|x — z| — > oo) one has: 

57^'"(x,y;z)~5o(x,2/), (9) 

where 5*0 denotes the free propagator. Typically, corrections to eq. (^) lead to small contributions and will be 
neglected in what follows. Once the propagator has been calculated, one can in principle evaluate any correlation 
function in the single-instanton background. 

Now, let's turn to the realistic vacuum of QCD. Here, any configuration with a non-zero net topological charge would 
be highly disfavored by the small value of the 0-angle. Therefore, one is lead to picture the vacuum as an ensemble 
with equal density of instanton and anti-instantons. If the vacuum is dilute enough, the classical back-ground field 
can be approximatively taken to be a superposition of separated instantons and anti-instantons [ p^ : 

A,{x, = ^ {nl },) + ^ A^^ix, },), (10) 

/ A 

where denotes the set of all collective coordinates. 

The propagator in such back-ground field can then be evaluated as follows ||l^. Let's consider the expansion: 

5 = 5o + ^o4^o + ^o45o45o + ..., (11) 

where integrations over the positions of each background field insertion is understood. The series ( |ll| ) can be re- 
arranged so that all terms depending on the collective coordinates of one instanton field only are summed up first, 
followed by all terms depending on two instantons and so on. One gets: 

S = 5o + ^(5/ - 5o) + Y,{Si - So)S,^{S,j - So) + (12) 
I i^j 

where Sj denotes the full propagator in the field of the instanton / so, in the approximation (^ one has: 

(S*/ - 5o)jj(x,?/) ~ : , (13) 



where we have dropped all collective coordinates indices. Inserting ( |13| ) in (12) and dropping also all spinor indices 
we get: 



I m 



E 



( I d^z^l,{z){^^. +^m)i^Jo{z) -imSjj) ^^^^ + (14) 



where —imSi^j has been added in order to relax the J ^ I constraint in the summation. All the terms, starting from 
the second on, form a geometrical progression, which can be re-summed to give: 



S{x,y) c S'{x,y) + J2^oiix) ( \ ) ^jiy), (15) 



I, J 

where T/j denotes the overlap matrix in zero-modes subspace 



Tij = / d''z^p\z)i{ip)4,iz)j. (16) 



In ( |15D , the zero- mode part the quark propagator is approximatively written as a bilinear form in the space spanned 
by the quark zero-mode wave functions. From (Q) it follows that the contribution coming from all the terms in the sum 
associated to instantons very far away from the points x and y will be negligible. In particular, the biggest term in 
( |l5| ) is associated to the closest instanton, /*. Such instanton is dominating if the average of the correlation function 
calculated retaining only the (/*j/*) term in ( p^ is much larger than the average of the same quantity calculated 
from all other terms in the sum (|l5[). Notice that this is a much weaker assumption than demanding 

V'o/»(x) ( \ ) ^J,.(2/)» J2 ^o/(x) ( \ ) ^ljiy\ (17) 
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for each configuration. 

Let us summarize the framework developed so far. First of all, the inverse matrix (y) jj contains all the information 
about the particular configuration of the instanton ensemble. In order to evaluate correlation functions, one needs to 
average over all possible configurations. Since contributions from distant instantons are suppressed by their zero- mode 
wave functions, one expects correlation functions with the highest number of zero-modes to be most influenced by the 
leading- instanton /* . If it is possible to retain only the contribution from /* , the global properties of the ensemble 
are present in the matrix element ■ 

As it was suggested long time ago by one of us ||l|, one can represent collective contribution of all other instantons, 
by introducing an effective mass associated to quark propagating in the zero-modes. In other words, one assumes that 
for |a; — y| < l/m, the quark propagator can be written as: 



S{x,y) 



^o(a;)'0j(y) 



(18) 



With such propagator all quark correlation functions in the instanton back-ground could be evaluated simply by 
computing all relevant Feynman diagrams and then averaging over the instanton collective coordinates [|9) . 
More specifically, in Random Instanton Liquid Model (RILM) one introduces a model instanton density n{p), 



ni{p) := ni d{p), 



(19) 



where nj = ua — \f'ni ^ and d[p) represents the instanton size distribution. The latter is schematically taken to be: 



d{p) = 5{p- p). 



(20) 



with p ~ 1/3 fm. This approach has the advantage to be considerably simple and was also proven to be quite 
phenomenologically successful [|[ |^. However, we show below that the effective mass defined in (|l^) is a quantity 
quite different from its naive estimate (Q). 

In order to clarify the statement, let us first consider the quark condensate: 

=< 0\Tru{x)uix)\0 >=< TrS{x,x) >, (21) 

where, in general, the average is done over all possible gauge field configurations. In the SIA is easily evaluated: 

-2p2 



<Q\u{x)u{x)\Q >= / d'^z / dpnd{p) 



[{z — x)^ + p'^Y ■"'^ "^ti 



(22) 



where, for reasons that will become clear shortly, we have denoted with niuu the quark effective mass and n := m + UA- 
After performing the integrations one finds: 



(23) 



for any normalized d{p). 

Now, repeating the same calculation in the full liquid 



gives: 



Xuu = { Tr 



I.J \^ J IJ 



(24) 



where, again, the average is made over all possible configurations of the ensemble. A comparison between (23) and 
(H) gives: 

(25) 



Tr 



Let us now consider another quark condensate: 

Xuudd :=< 0\Tr[u{x)u{x)]-Tr [d{x)d{x)]\0 >=< [TrS{x,x)]^ > 
Such condensate receives double contribution from zero-modes. In the SIA one obtains: 

Xuudd ^ I dp dip) ^ " , 



5 TT^ 



(26) 
(27) 



idd 
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TABLE L Quark condensates evaluated in the full instanton ensemble and from the leading-instanton, only. 



condensate 


complete calculation 


LI 


\uudd 


(-2S2±5MeVf 


(-198±lMeVf 


(310 ± 7MeV f 


(309 ± 3AieV f 



where we have now denoted with rriuudd the quark effective mass. 

Comparing, as before, with the result of full liquid calculations leads to: 



dpd{p) 



5 TT^ p'^ 



1 



(28) 



Now, if the effective mass is universal, (m„„)^ 



m 



Tr 



uudd^ it would imply: 

2 



5 TT^n 



I dpd{p)X 



5 Tr'^np'^ 



(29) 



where we have used the ansatz ( pO| ) [21 . Some comments on eq. ( |29| ) are in order. First of all, in general quark 
condensate is rather inhomogeneous, and for parametrically dilute instanton ensemble this ratio is small. However, 
with empirical diluteness it happens to be not so small, about 0.6 . In principle, by measuring the left-hand side and 
right-hand-side of (^9|) on the lattice separately, one can estimate the accuracy of the universality of the effective 
mass. 

However, since different configurations and even points have different leading instanton, the corresponding value 
T/./. fluctuates, and the average of its different powers in general leads to different effective masses. (This effect 
should not be confused with the inhomogeneity of the condensates discussed above.) Let us define a parameter i?,„, 
such that Rm = 1 means universal mass (m„„)^ = mf^^^j^: 



Rm :— 



Tr 



5TT^p^n(\Trj:jji;oi{x){^)jji;lj{x) 



(30) 



III. NUMERICAL STUDY OF THE SINGLE INSTANTON APPROXIMATION. 

In general, reliability of the SIA depends on the vacuum diluteness. In this section, we want to establish whether 
the QCD vacuum with realistic density is actually dilute enough for the leading-instanton to be dominant, at least 
for some observables. 

For this purpose we have performed numerical analysis of several correlation functions, measured in the random 
instanton liquid model. In such ensemble, the vacuum expectation values are obtained by averaging over configurations 
of randomly distributed instantons of size p — 1/5 fm. The contribution from the leading-instanton is evaluated by 
retaining only the largest term in (^5|), for each configuration. 

We begin by considering two quark condensates Xuu and Xuudd, introduced in (|2^) and (p6[). We will show later 
that they represent all generic observables which receive contribution from one and two zero- mode propagators, 
respectively. 

In this calculation we average 5000 configurations of 20 instantons in a box of volume 3.4 x 1.8'^/m'*. The results 
of this simulation are presented in table || 

From these results one can see how the accuracy of SIA (keeping only the closest instanton) depends on the particular 
matrix element being evaluated. Naturally, the accuracy increases with the dimension of the operator involved, because 
it diminish the contribution of distant instantons. Specifically, SIA for dimension-six local operators which receive 
contribution from two zero-mode propagators agree with full calculation within a few percent. On the other hand, 
prediction for operators/correlators with only one zero- mode propagators are not really accurate: the error in quark 
condesate is large ( > 35% ). 
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TABLE II: Universality parameter, Rm and the effective masses evaluated in the RILM and in the IILM 



Quantity 


RILM calculation 


IILM calculation 


Rm 


0.4 


0.2 


UOMeV 


177MeV 


\/^uudd 


65MeV 


QlMeV 



Next we consider two-point correlation functions. This allows us to determine the scale at which the closest instanton 
is no longer dominant. At this purpose we have measured the pion pseudo-scalar two point function, 

P(x) :=<0|J5(x)jt(0)|0>, (31) 

where, 

J^^x) := u{x)^5d{x). (32) 

This particular choice is motivated by the fact that such correlation function is known to receive maximal contri- 
bution from quark zero-modes One expects many instantons effects to become important for larger than the 
instanton size and smaller than the typical distance between two neighbor instantons: 

l/3/r7i < |a;| < 1/to (33) 

Results of simulations including the contribution from all instantons and form the leading-instanton only are reported 
in figure (|l|). One can see that the agreement is lost for rather large values of {\x\ ^ 0.6/m). 

In the last section we saw that the SIA does not only assume leading-instanton dominance, but involves some 
effective mass parameters, which collectively describe the effects of all other instantons (see eqs. ( p5|) and (p^). In 
the next section we shall determine numerically such parameters. 



IV. NUMERICAL STUDY OF THE QUARK EFFECTIVE MASS PARAMETERS 

In section || we argued that the universality of the effective mass, which collectively describes the effects of all 
non-leading- instantons, can be put in relation to the fluctuations of the quark condensates through eq. ([29|). 

Obviously, the accuracy of calculations in the SIA depends on the value of R„i (defined in (|3^)) in realistic ensembles. 

We have have evaluated Rm and the corresponding effective masses, numerically in the random instanton liquid 
and in the interacting liquid (for a review of these ensembles see |^). Our results are summarized in table ||; 

These results show that, in the instanton vacuum with realistic density, the universality does not hold 

(34) 

This implies that an effective mass extracted from the quark condesate can not be used in calculations involving more 
than one zero-mode propagator. 



On the other hand, the results of numerical simulations presented in section III have shown that matrix elements 
involving only one zero-mode propagator (like the quark condensate) can not be reliably evaluated in the SIA, simply 
because the leading instanton is not dominant. As a consequence, one is forced to consider only correlation functions 
involving at least two such propagators and therefore iriuu is of no practical usefulness. 

In more general terms, one may address the question whether the effective mass parameter depends on the particular 
correlation function being evaluated. If so, this feature would spoil much of the predictive power of the SIA. In such 
pessimistic scenario the SIA would only allow to work out the functional expressions of small-sized correlations, but 
not their overall normalization. However we will show that the effective mass parameters depend essentially on the 
number of zero-mode propagators involved, and that rii^^^^ is in a way universal for a number of applications. In this 
case, SIA is predictive including the normalization. To check that we have extracted from the analysis of several 
hadronic two-point functions evaluated in SIA and in the liquid. In particular, we considered the pion pseudo-scalar 
the scalar diquark and the a nucleon scalar correlation functions: 

P{x) = < 0|J5(a;)J5t(0)|0>, (35) 
D{x) = <0\J^,{x)J^lmO>, (36) 
N{x) = <0|rr[r;(a:) 77(0)74] |0>, (37) 
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TABLE III: Estimates of the quark effective mass from several correlation functions. 



Correlation function 


ml [MeV^] (RILM) 


ml[MeV^] (IILM) 


Xuudd condensate 
pion pseudo-scalar 
diquark scalar 
nucleon scalar 


(65)2 


(91)2 


(65)2 


(105)2 


(69)2 


(105)2 


(67)2 


(105)2 



TABLE IV: Determinantal masses evaluated in the RILM and in the IILM as compared to mi and m^, defined in section 



mass 


RILM calculation 


IILM calculation 


mi 
mdet 


1207^6^ 


177MeV 


63MeV 


W2MeV 


ml 
■maet 


(65MeV)2 


(103MeV)2 


(64Mel/)2 


(103MeV)2 



where, 



r]a{x) 



u(x)-f^d{x). 



^ah c 



Ub{x) C75 dc{x). 
\u''{x)C-i,u\x))ul{x). 



(38) 
(39) 
(40) 



All these correlations function are known to receive contribution from two propagators in the zero-mode. 

The comparison between results obtained in the SIA, in the random instanton liquid model (RILM) and in the 
interacting instanton liquid model (IILM) are reported in figs. (^), (||) and (jj). The corresponding values for 



are presented in table (HI). These values are indeed rather different from the traditionally adopted estimate 



TO* = llQMeV, extracted from the quark condensate. 

The general reason why these masses are rather small is the following. Instantons have fluctuating strength of 
interaction with others in the ensemble: some of them are "hermits" and have small matrix elements in the corre- 
sponding entries of the overlap matrix T . As in all expressions we average the inverse of this matrix, the contribution 
of such "hermits" is enhanced. This lowers the value of the effective masses. Furthermore, because random ensemble 
of RILM has more such "hermits" , as compared to IILM (where the fermionic determinant in the statistical weight 
suppresses them), these masses are smaller in RILM as compared to IILM. Such discrepancy reflects the fact that the 
two ensembles give actually quite different correlation functions . 

From these results we conclude that ra\ seems to be a universal parameter, describing the collective many-instanton 
effects. 

It is important to know what value of to| is suggested by the available phenomenology. As before, we chose to 
consider the pion pseudo-scalar correlator, because it receives maximal contribution from instanton zero-modes. The 
traditional "pole-plus-continuum" model for the spectral decomposition of P{x), gives [llF 



|13 



3 /"^ 

P{x) = XlD{m^;x) + — - / ds sD{y/s; x), 



(41) 



where D{m\ x) is the scalar propagator, so is the threshold for the continuum ( ^Jsq ~ l.QGeV) and the pseudo-scalar 
decay constant A^r is given by: 



f Tfl 

rriu + nid 



(42) 



We determined m^, by fitting the SIA prediction to the phenomenological curve obtained from (pl|). We found (see 
fig- |): 

mlj,hen. = {SQMeVf. (43) 

To further check the approach, we have evaluated the scalar proton two-point function N{x), using the value ( |4^ ) 
and we have compared with the phenomenological curve (see fig. (^) [p3[). In summary: with this value we obtained 
good very good agreement with phenomenology and therefore we suggest that (^3|) should be used f the applications 
of the SIA, when two zero-mode propagators are involved. 
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V. EVALUATION OF AN EFFECTIVE MASS IN THE FERMIONIC DETERMINANT 



The propagator is not the only place where the Dirac operator appears: the QCD statistical sum contains its 
determinant, appearing in power given by the number of light quark flavors Nf. If one considers the academic 
vacuum with only one instanton, this determinat contains the product of "current" quark masses for all quarks Q . If 
this would be the final answer for the instanton density, instanton effect would be strongly suppressed by their small 
values. 

However, in physical vacuum there are sufficiently many instantons to break chiral symmetry and produce non-zero 
quark condesates and effective quark masses, which substitute for much smaller "current" masses and make instanton 
effects significantly stronger. The interplay between these effective masses and current quark masses is especially 
interesting for strange quark, since the former and the latter, mg, are of comparable magnitude. This issue has been 
studied e.g. in recent paper [l6| , where it was concluded that the usual additive formula for total effective quark mass 
of the strange quark M*°* = nief / {mg = 0) + is wrong, and the true value of Af*°* is not very different from that 
for u,d quarks because me//(TOs) strongly decreases with m^. 

Apart of the role of strange quark mass in general, there is also a general issue of correct connection of units and 
vacuum parameters (with the instanton density being one of them) for QCD with different number of flavors. (For 
example, between no-quark or quenched QCD and the physical world.) In order to study all of this, it important to 
know what is the absolute magnitude of the fermionic determinants in the instanton-based vacuum models considered. 
Some of those are repoted in this section. 

In the instanton-based model context, the fermionic determinant is usually represented by the determinant of the 
overlap matrix T (see description e.g. in |^) in the zero mode subspace. After averaging over appropriate ensemble, 
one can define the so-called "determinantal masses" : 

<,:^^<«^, .^1,V.^ (44, 

< > 

where index i refers to number of flavors and N denotes the number of instantons. Their values tell us how much the 
presence of fermions reduces the instanton density, compared to the same ensemble without them. 

Originally, in an estimate for the determinant effective mass was estracted from the averaging of the 't 

Hooft Lagrangian assuming factorization of quark condesates, and using the same m* — 170 MeV. If so, each flavor 
reduces instanton density by the factor m*p w 0.28. As we will see shortly, the corresponding reduction factor is 
actually even smaller. 

In principle, there is no reason why the values of m^iet and rn^^^ should agree with mi and mj, deflned in the 
previous section: we now averahe positive rather than negative powers of the overlap matrix. 



We have evaluated the determinantal masses in the RILM and in the IILM. Results are reported in table [V. Some 
comments are in order. First of all note that, in both ensembles, the values of m^g^ turn out to be quite consistent 
with the values of Furthermore, the fluctuations of the determinantal mass, mdet and [mdet)"^ are very small: 

■m^det ~ {"^detf ^ ml - [mif , (45) 

implying essentially that mi is inconsistent with mdet- This fact could have two possible explanations. On the one 
hand, one could argue that mi is a somewhat ill-defined parameter, because the SIA can not be used to evaluate 
quark condensate. On the other hand, one could observe that larger fluctuations for the effective masses defined in 



IV should not be surprising, since such parameters appear always in denominators of SIA calculations. 



VI. CONCLUSIONS AND OUTLOOK 



Summarizing our study of the SIA approximation in QCD, we first notice that this approach has been related to 
the theory of the full ensemble and all the effective parameters previously loosely called "effective masses" are defined. 
All of them describe different aspects of collective interaction between the "leading" instanton (the closest to the 
observation points) and all others, and related to the overlap matrix T . Different effective mass values simply follow 
from different ensemble averaging. In particular, the factor appearing in SIA calculations with one propagator 
in the zero-mode, does not correspond to the square root of the factor appearing when two such propagators are 
involved. 

We have made numerical simulations in the RILM and IILM and found that the contribution of the leading-instanton 
actually dominates all condensates of operators of dimension six or more, as well as short-distance correlation functions 
(|a;| ^ 0.6/m). This however is true only for correlation functions with at least two zero-mode propagators involved. 
Earlier estimates extracted from the quark condensate are not accurate. 
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Furthermore, the parameter -\ is approximatively universal for several correlation functions with two zero-mode 

7712 

propagators involved. We have also extracted a phenomenological estimate of its value from the analysis of the pion 
pseudo-scalar correlator. We found m^pi^^^ — {86MeV)'^ , much smaller than the value originally obtained from the 
quark condensate. Our new value should be used in many applications of the SI A. 

Finally, we have compared our estimates for the effective mass parameters mi and m^, with the measurements of 
the "determinantal" masses, introduced in We observed substantial agreement between m| and m^^j both in the 
RILM and the IILM, but diffe 3 implies that light quarks 

are about twice more effective 




FIG. 1: Pion pseudo-scalar correlation function in the RILM, normalized to the same correlation function in the free theory. 
The solid line corresponds to th' ton contribution. 
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FIG. 2: Pion pseudo-scalar two-point function normalized to the same correlation function in the free theory. The open circles 
(squares) represent RILM (IILM) points, the dashed lines represent SIA calculations with masses given in table |l| and the 
dotted line is the phenomenological curve obtained from the spectral decomposition. 



10 



e 2^ 



B —-e— 



0.2 



0.4 

x,[fm] 






RILM 


□ 


IILM 




SIA 



0.6 



0.8 



FIG. 3: Diquark scalar two-point function normalized to the same correlation function in the free theory. The open circles 
(squares) represent RILM (IILM) points and the dashed lines represent SIA calculations with the effective masses given in table 
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FIG. 4: Nucleon scalar two-point function normalized to the same correlation function in the free theory. The open circles 
(squares) represent RILM (IILM) points and the dashed lines represent SIA calculations with the effective masses given in table 
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For simplicity, we shall often use the term "instanton" to denote instantons and/or anti-instantons. 
For sake of simplicity, we are explicitly dropping all collective coordinates, except the instanton position Zi; moreover the 
use of the singular gauge is assumed everywhere. 

Notice that for all gauge invariant matrix elements, the average over the color orientation is trivial. 
Here, we have neglected all small current quark mass terms in 1/T. 

Alternatively, we repeated the calculation using a parameterization of the lattice measurements of d{p), which is peaked 
about somewhat higher values of p (p ~ 3.9). Both calculation give basically the same result. 
Here we have averaged on 5000 configuration of 256 instantons in a 4'*/m'* box. 

Notice that, in this case, the pole contribution depends on the nucleon "decay-constant" As, which is not known experi- 
mentally. We have therefore used an estimate (A^ — 2.5fm~^) which reasonably agrees with several model calculations P] 



